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1. Introduction 

Supersymmetry plays a central role in theoretical models for elementary particle physics be- 
yond the Standard Model. Therefore it is important to gain knowledge about the properties of 
supersymmetric theories. Much of what is known about supersymmetric models is based on tree- 
level considerations or comes from perturbation theory. However, various important characteristics, 
like the existence of mass-degenerate supermultiplets of particles, are of a non-perturbative nature. 
Therefore it is desirable to study them by means of non-perturbative methods. 

The simplest supersymmetric model including gauge fields is the supersymmetric Yang-Mills 
theory (SYM). It describes interacting gluons and their supersymmetric partners, the gluinos, which 
are Majorana fermions in the adjoint representation of the gauge group S\J(N C ). The (on-shell) 
Lagrangian of SYM in Minkowski space is 



Jz^ = Tr 



X -F^ v + l -Xy^D,X-^XX 



(1.1) 



where X is the gluino field, the gluon field, F^ v the non-Abelian field strength, and 

D li X a = d li X a + gf ab cA b ll X e (1.2) 

denotes the gauge covariant derivative in the adjoint representation. The gluino mass term breaks 
supersymmetry softly. 

SYM is similar to QCD in some respect [|l|]. It is asymptotically free and is assumed to show 
confinement. The "physical" particles are bound states of gluons and gluinos, and if supersymme- 
try is unbroken, they would form supermultiplets. The non-perturbative properties one would like 
to investigate with the lattice simulations include: (1) the spontaneous breaking of chiral symmetry, 
Zin c —> Z2, that manifests itself in the non- vanishing vacuum expectation value (XX) / 0, (2) the 
confinement of static quarks, indicated by a linear rise in the static quark potential, which is an evi- 
dence for the confining nature of the theory, (3) the spectrum of bound states, that can be compared 
to the predictions of low energy effective actions, (4) whether supersymmetry is spontaneously 
broken, (5) the restoration of SUSY in the continuum limit of a lattice regularisation. 

The particle content of SYM is expected to consist of colour neutral bound states of gluons 
and gluinos, forming supermultiplets. Based on effective Lagrangeans, it has been predicted [Q, ||] 
that the low-lying particles form two chiral supermultiplets, each consisting of a scalar, a pseu- 
doscalar, and a fermionic spin 1/2 particle. One of them contains a 0~gluinoball (a-T}' ~ Xy$X), 
a + gluinobail (a-/o ~ XX), and a spin 1/2 gluino-glueball (% ~ a^ v Tr(F^ v A)), the other one 
a 0~ glueball, a + glueball, and a gluino-glueball. Both supermultiplets contain an exotic particle 
state called gluino-glueball, which is a spin 1/2 Majorana fermion. Such a bound state containing 
a single fermion does not occur in QCD, but analogous particles exist in models similar to QCD 
with an arbitrary number of quark flavours in the adjoint representation. 

We study SYM non-perturbatively in the framework of regularisation on a lattice. Supersym- 
metry is generically broken in any non-trivial theory on the lattice [B]. The particle spectrum and 
the supersymmetric Ward identities can show how it is restored in the continuum limit or whether 
there is a possible remnant breaking of SUSY by non-perturbative effects. The existence of super- 
multiplets is an important signal for the supersymmetric limit of the theory. 
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Previous work by our collaboration on SYM on the lattice, concerning the nonperturbative 
items mentioned above, is reported in [f| and references given there. The results obtained there 
have not yet shown the expected degeneracy of the fermionic and bosonic masses. The mass of the 
gluino-glueball appeared to be larger than the other masses of its lightest possible superpartners. 
However, the masses were obtained at a fixed lattice spacing and without a detailed analysis of 
the finite size effects. Below we shall discuss our recent calculations, which indicate that the 
influence of the finite lattice spacing is larger than expected and provides a possible source of the 
supersymmetry breaking in the simulation. 



2. SUSY on the Lattice 

In our work we employ the formulation of SYM on a lattice by Curci and Veneziano [^], 
where the gluinos are described by Wilson fermions in the adjoint representation. SYM has also 



been investigated with domain wall fermions [|8j, g, |10|] and overlap fermions flllQ, which, however, 
require significantly more computing resources than Wilson fermions for large lattice volumes and 
small lattice spacings. 

The lattice action of SYM in our setup is 

Sl = j3£ (l - ^"Re TrU^ + 1 £ A,(D)^ , (2.1) 
where D is the Wilson-Dirac operator 

(P)x,a,a;y,bfi = 5xy5a,b8a,p 
4 

L [C 1 _ faW ( V M M WWv + 0- + Yn)a,P (V£ (* ~ M) WSt-w] , (2.2) 

and (y^{x)) a b are the gauge link variables in the adjoint representation. We are currently consider- 
ing the gauge group SU(2) with generators T a = X a /2, in which case 

%(x)U = 2Tr (U^T a U m T b ). (2.3) 

The hopping parameter K is related to the bare gluino mass via K = 1/ (2m g + 8). 

In order to reduce lattice artifacts, in our simulations we actually use the tree-level Symanzik 
improved gauge action, and one level of stout smearing applied to the link fields in the Wilson- 
Dirac operator. 

From considerations of chiral and SUSY Ward identities it is expected [[7|] that a fine-tuning of 
the bare gluino mass parameter (i.e. k) in the continuum limit is sufficient to approach the chiral 
symmetry and supersymmetry of the continuum theory. This tuning is most efficiently done by 
means of the mass of the unphysical adjoint pion (a-7i). This particle is the pion in the correspond- 
ing theory with two Majorana fermions in the adjoint representation. The correlator of this particle 
is the connected contribution of the a-T]' correlator. The a-7T is not a physical particle in SYM. 
However, it can be defined in a partially quenched setup, in the same way as for one-flavour QCD 



[12]. On the basis of arguments involving the OZI-approximation of SYM M\, the adjoint pion 
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mass is expected to vanish for a massless gluino. The corresponding value of K c is most easily 
obtained from the dependence of the a-7T-mass on fc. 

The numerical simulations are done using a polynomial hybrid Monte Carlo (PHMC) algo- 
rithm with a two step polynomial approximation and reweighting fljjj §]. The functional integral 
over Majorana fermions yields the Pfaffian Pf(CD) = sign(Pf (CD)) %J det(CD) of the antisymmet- 
ric fermion matrix multiplied by the charge conjugation matrix C. The PHMC implements the 
positive weight y / det(D), while the sign sign(Pf(CD)) is taken into account by reweighting. The 
reweighting factors are obtained from the number of negative real eigenvalues of D. 



3. The particle spectrum on the lattice 

The most interesting particles are the possible candidates for the lowest lying multiplets. These 
multiplets consist of a fermionic, a scalar, and a pseudoscalar particle. On the lattice a mix- 
ing between possible operators representing the corresponding quantum numbers is expected. As 
representations of the bosonic scalar and pseudoscalar particles we have measured glueballs and 
gluinoballs. The fermionic particle is represented by the gluino-glueball. 

The correlation function of the ++ glueball is calculated by smeared plaquette-plaquette cor- 
relations. Variational smearing techniques (APE and HYP smearing) are used to get a reasonable 
signal for the particle mass. Unfortunately the signal for the v glueball turned out to be insuffi- 
cient for a determination of the mass with our current statistics. 

The gluinoballs are similar to flavour singlet mesonic operators in QCD. Hence they contain 
disconnected contributions and require the computation of all-to-all propagators. In a graphical rep- 
resentation of the correlators the disconnected contributions show up as two disconnected fermion 
loops, e. g. for the a-rj' 



{l(x)y,X{x)X(y)y 5 X(y)) = {\J\J-2*Q* >• (3-D 

From the computation of flavour singlet mesons in QCD it is well known that the all-to-all propa- 
gators introduce additional statistical noise in the computation of an observable. We observe this 
effect in our simulations when we compare the larger statistical errors of the a-Tj' with the small 
error of the a-7T. We have implemented and tested several techniques to determine the all-to-all 



propagators H14Q. In our recent simulations a combination of the truncated eigenmode approxima- 
tion and the stochastic estimator method has been used. Our truncated eigenmode approximation 
is implemented in such a way that from the eigenvalues reweighting factors can be obtained to im- 
prove the polynomial approximation of the PHMC algorithm. The stochastic estimator method was 



improved with iterations of a truncated solver, see [ |15[ ] for the details of these methods. With these 
techniques the a-Tj' mass can be obtained with a reasonable precision, while the determination of 
the a-/o mass is still not satisfactory. 

In contrast to the difficulties in the bosonic sector, the mass of the fermionic gluino-glueball 
can been obtained most reliably. The corresponding lattice operator is a' IV Tr[F uv A], where the 
field strength on the lattice is represented by the clover plaquette. A combination of APE (applied 
on F^ v ) and Jacobi smearing (applied on A) leads to a further reduction of the noise for this signal. 
Plots of the effective masses for this particles can be found in Fig. [l[ 
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Figure 1: Examples for the effective masses of the particles measured in supersymmetric Yang-Mills theory. 
Time slice distances t,t + a are used to compute the effective mass m e g at distance t assuming a single cosh 
(sinh) behavior. The mass and t are given in units of the lattice spacing a. The lattice size is 24 3 x 48. The 
parameters of the simulations are j3 = 1.75 and K = 0.1490/0.1492. 

4. Finite size effects and lattice artifacts 



In previous work of our collaboration a rather large gap between fermionic and bosonic masses 
has been observed [J3p - These simulations were done at j8 = 1.6 and lattice sizes up to 24 3 x 48. 
It would, however, be premature to draw conclusions from this about possible supersymmetry 
breakings in continuum SYM. There are several limits involved in the determination of physical 
results for this model, namely the infinite volume limit, the continuum limit, and the chiral limit, 
which have to be taken in this order. With present resources it is not possible to carry out this 
program. In recent simulations we have, however, been able to obtain some estimates for the 
relevance and influence of these limits. 

We have performed a detailed investigation of the finite size effects in simulations of 8 3 x 16, 
12 3 x 24, 16 3 x 36, 20 3 x 40, 24 3 x 48, and 32 3 x 64 lattices at /3 = 1.75 [g6||. For the asymptotic 
dependence of masses on the box size L we have taken [[n|, [TsJ] 



m(L) zmiQ + CL 1 exp (—amoL) 



(4.1) 
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Figure 2: The influence of the finite box size L on the mass gap between the bosonic a-rj' and the fermionic 
gluino-glueball. These results were obtained at a fixed value of the bare parameter K = 0.1490 (j3 = 1.75). 
The mass is, as indicated, shown in units of r* 1 , and L in units of ro/0.5, where ro is the Sommer parameter. 
When ro is set to its QCD value, the units of L would correspond to fm. 

for large values of L. Results have been obtained at fixed bare parameter fc and at a fixed mass of 
a-7T. As explained above, the best signal is obtained for the masses of the gluino-glueball and the 
a-T]'. The difference between their masses can be taken as an estimate of the mass gap between 
the particles of a supermultiplet. The dependence of the mass gap on the finite box size is shown 
in Fig. g. 

The mass gap is clearly increased by the influence of the finite volume. The supersymmetry 
breaking terms, that are present due to the lattice discretisation are hence increased by the finite 
size of the lattice. However, this effect becomes rather small already at moderately large lattice 
sizes. Therefore, finite size effects cannot be the source of the observed SUSY-breaking effects. 

We have also been able to estimate the influence of the finite lattice spacing from simulations 
at a second larger value of j8. The difference between the results at j8 = 1.6 (a 0.08 fm) and 
/3 = 1.75 (a « 0.06 fm) are shown in Fig. [|. They indicate that indeed the lattice artifacts are the 
most relevant source for a bias in the spectrum of the particles. 

Possible alternative ways to reduce the lattice artifacts are currently being explored. For ex- 
ample, an increased smearing in the fermionic part of the action seems to reduce their influence. 

5. Recent results 

The results of our investigations of finite size effect and lattice artifacts are the basis for our 
current simulations. The influence of the finite volume seems to be small. A small lattice spacing 
is important to keep the influence of the lattice artifacts under control. This result, obtained from 
the gluino-glue and the a-T]' mass, is already quite convincing. For a complete investigations of 
the effects it is, however, indispensable to get information also form the other particle states. To 
distinguish the small mass gap at /3 = 1.75 from the statistical and systematic errors, a good signal 
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Figure 3: The mass of the a-T}' and the gluino-glueball at two different j3, where the larger value corresponds 
to the smaller lattice spacing. The masses are shown as a function of the squared mass of the adjoint pion 
(mji). The results from the two different lattice sizes 24 3 x 48 and 32 3 x 64 are found to be consistent. All 
masses are given in units of the Sommer parameter tq. 
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Figure 4: The masses of all particles at /3 = 1 .75 (a ~ 0.06 fm) obtained on a 32 3 x 64 lattice (L= 1.8 fm) 
with a statistics of about 5000 configurations at each point. All masses are given in units of ro. 



is needed for all observables. Before we test further improvements of the action to reduce the lattice 
artifacts we have, therefore, made improvements in the measurement of the glueballs and the scalar 
mesons. Especially for the glueballs most important for an improvement is to increase the statistics 
significantly. Our first investigations were done at a lattice size of 32 3 x 64 (L = 1.8 fm), see Fig. |j. 
In our investigations of the finite size effects we found that these results are consistent with those 
at a lattice size of 24 3 x 48 (L = 1.35 fm ). Based on these findings we have chosen the smaller 
lattice size and increased the statistics from about 5000 configurations at each K to around 10000. 
A preliminary summary of these data is shown in Fig. ||. 
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Figure 5: The masses of all particles at /3 = 1 .75 (a « 0.06 fm) obtained on a 24 3 x 48 lattice (L = 1 .35 fm) 
with a statistics of about 10,000 configurations at each point. All masses are given in units of tq. 



6. Conclusions 

We have investigated the influence of the various simulation parameters on the unexpectedly 
observed mass gap between the particle states of a supersymmetric multiplet. 

The finite size effects increase the mass gap as found from the splitting of the a-Tj' and gluino- 
glueball mass. The supersymmetry breaking terms that are introduced by the lattice discretisation 
seem to get finite size collections that enlarge their influence. However, this effect is found to be 
small compared to the statistical error at moderate lattice sizes above around 1.2 fm. 

The finite lattice spacing was the most relevant source of supersymmetry breaking in the previ- 
ous results of our collaboration. The mass gap at our current, smaller, lattice spacing is significantly 
reduced compared to the previous results. In view of the statistical errors of the observables, it is 
already difficult to obtain the mass gap. 

To determine the remaining mass gap with enough precision and to answer the question, 
whether it persists in the continuum limit, it is necessary to reduce the statistical errors of the 
observables. One important parameter for such an improvement is the statistics. An investigation 
of the improvements of the observables is the main task of our current simulations. 

Improvements of the action, like an increased level of stout smearing or the inclusion of a 
clover term might also help to get a more precise answer about the relevance of the mass gap in the 
continuum limit. We have started to investigate these improvements. 
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